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Abstract 



Arrhenius-type plots for multifragmentation process, defined as the trans- 
verse energy dependence of the single-fragment emission-probability, (ln(l/pb) 
vs 1/y/Et), have been studied by examining the relationship of the parameters 
Pb and Et to the intermediate-mass fragment multiplicity (n). The linearity 
of these plots reflects the correlation of the fragment multiplicity with the 
transverse energy. These plots may not provide thermal scaling information 
about fragment production as previously suggested. 
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About a hundred years ago, the Swedish chemist Svante Arrhenius discovered that the 
rate of chemical reactions increases with temperature [1]. Specifically, the chemical reaction 
rate constant (k) is related to the absolute temperature (T) by 

k oc exp{-E a /T) (1) 

where E a is the activation energy of the chemical reaction. The linear relationship between 
ln(k) and 1/T is widely known as Arrhenius plot in chemistry and is associated with thermal 
equilibrium processes. 

Recently, Arrhenius-type plots have been used to study the statistical [2-5] and dynami- 
cal [6] properties of fragment emissions in heavy ion reactions. In the intermediate incident 
energy range, between few tens of MeV to 100 MeV per nucleon, the production of interme- 
diate mass fragments (IMF, 3 < Z < 20), also known as multifragmentation, is an important 
decay mode of highly excited nuclear systems [7]. Calculations suggest that the fragments 
are produced in the phase co-existence region. Thus, understanding the mechanisms of frag- 
ment formation may provide an insight to the liquid gas phase transition in nuclear matter 
[7]- 

Several experiments [2-4] gave evidence that in collisions characterized by a given value 
of the total transverse energy of detected charged particles, 

£ t = 5> jS in 2 ^ (2) 

% 

the IMF multiplicity distribution may be fitted by a binomial distribution, 

TO ' 

P n = "17 " VyPbi 1 ~ (3) 

n\(m — n)\ 

where n is the IMF multiplicity and the parameter m is interpreted as the number of times 
the system tries to emit a fragment. The probability of emitting fragments can be reduced 
to a single particle emission probability pb for true binomial distributions; the binomial 
parameters m and pb are related to the mean and variance of the fragment multiplicity 
distributions according to: 

2 



(n) = mp b 



(4) 



ol = (n)(l-p b ) 
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The past investigations found a simple linear relationship between ln(l/pb) and 1/ \[E~ t (nu- 

clear Arrhenius-type plots), for several projectile-target combinations and incident energies 
[2-4]. By assuming a linear relationship between \[E~ t and temperature T and from the lin- 
earity of the observed ln(l/p&) vs l/y/E~ t plot, it has been inferred that a thermal scaling of 
the multifragment processes might be a general property [2-4]. In this picture, the observed 
"linear" dependence of ln(l/p) upon 1/T would be reflecting the 



dependence of fragment emission probabilities upon a common fragment emission barrier B. 

However, unlike chemical reactions, p and T were not measured directly in Refs. [2- 
4]. There, the validity of the Arrhenius-type plots relies on two assumptions: 1) that E t 
is proportional to the excitation energy E*, and therefore, should be proportional to T 2 , 
and 2) that pb obtained by fitting the fragment multiplicity distributions is the elementary 
emission probability p. This article will examine the above assumptions and investigate the 
underlying reasons for the linearity exhibited by the Arrhenius-type plots obtained in many 
systems. 

The assumption that temperature is proportional to the square root of the excitation en- 
ergy, T oc y/~E*, is valid for compound nuclei formed at low to moderate temperature. Some 
relationship between the transverse energy and excitation energy and therefore temperature 
may also be obtained for compound nuclei, provided adjustment is made for the Coulomb 
barrier and for neutron emission [8]. For the intermediate-energy heavy ion reactions such 
as 36 Ar + 197 Au at E/A = 35 to 110 MeV, where linearity of the Arrhenius-type plots 
have been observed, however, the final states contain fast particles emitted from the overlap 
region of projectile and target as well as delayed emission from projectile- and target-like 
residues. In particular, particle production from the overlap participant region dominates 



p oc exp(-B/T) 



(6) 
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in central collisions [9-12], and the transverse energy from this region is strongly affected 
by the collective motion [13]. There has never been any unambiguous experimental evi- 
dence supporting that \fE~ t is proportional to the temperature and in fact the evidence is to 
the contrary. Recent temperature measurements using both the excited states populations 
and isotope yield ratios show that the temperature dependence on the impact parameter 
determined from charge particle multiplicities is very weak [14-16], less than 1 MeV from 
peripheral to central collisions. Similar trends have also been observed for the Ar + Au 
reactions at E/A = 35MeV [17]. Since E t is strongly impact parameter dependent [18], 
these temperature measurements thus imply that E t is independent of temperature and the 
assumption that T oc \[E~ t is not valid. Based on this argument alone, the ln(l/pb) vs 1/ \fE~ t 
plot is not the Arrhenius plot analogous to that observed in chemical reactions. The break- 
down of the T oc y/E~ t assumption suggests that previous interpretation of thermal scaling 
on emission probabilities, charge distributions and azimuthal correlations [2-4, 20-23] should 
be re-examined. 

Next, we will examine the assumptions used to extract the fragment emission probability. 
There is no apriori reason for the emitted fragments to prefer binomial statistics or Poissonian 
statistics. In Poissonian statistics, the probability of emitting n fragments is 

P p (n) = -e"* (7) 

where A = (n) is the mean. The major difference between the binomial and Poisson dis- 
tributions is the ratio of the variance to the mean, a 2 / (n) where a 2 j (n) = 1 for Poisson 
distribution and < 1 for binomial distribution. It has been demonstrated that constraints 
from conservation laws reduce the width of the Poisson distributions to much less than 1 
[19]. For example, if charge conservation constraint is applied to a Poissonian distribution, 
Eq. (7) is modified to (see Appendix) 

P m (n,a) = ^e-V^- A ) 2 (8) 

where a is the charge constraint factor. 
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For small a, the mean fragment multiplicity for Eqs. (7) and (8) are nearly the same; 
(n) Rs A. Fig. 1 shows three modified Poisson distributions of Eq. (8) (solid and open 
points) for A = 3, 6 and 10 and a = 0.1. To illustrate that distributions such as Eq. (8) 
whose values of a 2 / (n) are less than 1 can be described by binomial distributions, we used 
Eqs. (4) and (5) to determine the binomial parameters m and pb whose values are listed in 
Fig. 1. The solid and dashed lines are binomial distributions of Eq. (3). The agreement 
between the two distributions is very good. However, in this context, m and pb are mainly 
fit parameters used to describe the modified Poisson distributions of Eq. (8) and pb is not 
an elementary emission probability. If the small values of a 2 / (n) reflect the constraints of 
conservation laws of Eq. (8) observed in Ref. [19], the reducibility of fragments emission to 
binomial distributions shown by Refs. [2-4] does not imply any fundamental significance for 
the parameters m and pb thereby extracted. 

Even though the ln(l/p fe ) vs \j\fE t plots constructed in heavy ion reactions are not true 
Arrhenius plots, analyses of many systems [2-4, 25] suggest that the approximate linearity 
of such plots may be universal. To explain this appealing systematics, we examine the 
correlations between the observable E t , parameter pb, and the fragment multiplicity n [25]. 
Since E t is obtained from the measured energies of both the light particles and fragments, 
Eq. (2), the energy and multiplicities are related by 

E t ^(N c -(n))E t LP +(n)El MF , (9) 

where Nc is the total charge particle multiplicities, E^ p and E IMF are the average transverse 
energy of a light particle and an IMF, respectively. Experimentally, the dependence of N c 
on (n) can be approximated by [9-12] 

N c = a' + b'{n) (10) 

where a' corresponds to the typical number of light charge particles emitted before any 
IMF is emitted and b' is the number of light charge particle emitted for each IMF emitted. 
There is some non-linear dependence of (n) at very high Nc, where (n) saturates for central 
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collision, and at very low Nc, where fluctuations in small value of n prevent a sharp cutoff 
in N c . Except for very large and small values of N c , Eq. (9) can be then rewritten into: 

E t = a + b(n) = b(a/b+ (n» (11) 

where a = a'E^ p is the threshold transverse energy associated with light particles emitted 
before any IMF and b = (b' - 1)E^ P + E IMF . 

The binomial fit parameter m is nearly constant as a function of the transverse energy 
E t [2-4]. Thus the plots of ln(l/p fe ) vs \/\fE t can now be reduced to ln(l/(n))[26] vs 
1/ \J a/b + (n) according to Eqs. (4) and (11). Fig. 2 shows the dependence of l/{n) and 
l/y/a/b+(n) for (n) ranging from 0.25 to 5.0, typical values observed in multifragmentation 
of heavy ion reactions [2-4,9-12]. For a/b = 0, the curve is concave and for a/b = 1, the 
curve becomes slightly convex. In the middle region, where a/b ~ 0.5, the curve is nearly 
linear. Thus the linearity of the Arrhenius-type plots merely reflects the correlation of the 
fragment multiplicity with itself when the value of a/b is about 0.5 which naturally arises 
from the intrinsic linear dependence of (n) on E t . 

To illustrate the self-correlation effect in nuclear Arrhenius-type plots, the published 
data for the Ar + Au collisions at E/A — 110 MeV [2, 4] are plotted as solid points in the 
left panel of Fig. 3. The solid line is the self-correlation of l/(n) and 1/^Ja/b + (n), scaled 
according to Eqs. (4) and (11) using the experimental determined values of m — 12 [2, 4], 
a/b = 0.5, b = 213 MeV [2, 4, 24]. The good agreement between the data and the self 
correlation confirms that the linearity observed in the Arrhenius-type plot mainly comes 
from the linear dependence of E t on (n) with a non-zero offset in Eq. (11). 

The non-zero value of a arises from Eq. (10) because a' is not zero. One would expect that 
a relation between 1 / (n) and 1 / \/W c similar to those shown in Fig. 2 should be observed. 
The right panel of Figure 3 shows the plot of l/(n) vs l/\/W c for Ar + Au reaction at 
E/A = 110 MeV [2,4,24]. The linearity is of the plot is comparable to most Arrhenius- 
type plots [2-4, 25]. Since (n) and Nc are much less affected by the energy resolution of 
the detection device, they are better observables than p b and E t used in the Arrhenius- 
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type plots. Figure 3 suggests that the Arrhenius-type plots contain essentially the same 
information as the much simpler plots of the IMF multiplicity ((n)) as a function of charge 
particle multiplicity (Nc) published in many studies [9-12]. 

In summary, recent temperature measurements show that temperature is nearly inde- 
pendent of impact parameter and therefore the transverse energy E t is not related to tem- 
perature. If y/E~ t is not proportional to the temperature, ln(l/p(,) vs l/V^t are n °t true 
Arrhenius plots. Without invoking the interpretation of fragment emission probability in 
binomial distributions or the temperature dependence of E t , the linearity of the Arrhenius- 
type plots can be reproduced from the linear dependence of E t on fragment multiplicity (n) . 
These plots carry the same information as N c vs (n) plots and do not provide additional 
information concerning the thermal scaling of fragment emissions. 

Appendix 

When many independent processes contribute to the multiplicity of certain particles, each 
with a small probability, then the multiplicity is expected to obey a Poisson distribution. 
Within the grand canonical ensemble, a Poisson distribution also follows. In general, the 
Poisson distribution is associated with a large overall system compared to the subsystem 
studied. However, in heavy ion reactions, the systems studied are always finite, constrained 
by conservation laws such as the overall energy, charge, and mass conservation within the 
participant region. 

To investigate the minimal effect of the above mentioned constraints on the multiplicity 
distribution of IMF, we first consider a situation where probabilities of emitting various 
individual particles are independent, with multiplicities being governed by a Poisson distri- 
bution in the absence of any constraint. When multiplicities of individual IMFs are governed 
by a Poisson distribution, then the overall multiplicity of IMFs is also governed by a Poisson 
distribution [27]. For simplicity, we next impose only a single constraint, that of the charge 
conservation, and examine changes in the multiplicity distribution of IMF. The constraint 
modifies the Poisson distribution to 
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Pm(n) oc £ ^ n (P P T M) E r oth (Z oth ) 5 z , Zoth+E p . (12) 

{n„:i/£lMF} r Z oth 

where is a Poisson distribution for fragment r, is the charge distribution of all 
particles other than IMF, and Z is the charge of the emitting source, Z = {Z ot h) +^ 1/ (n u )z u . 
If the system emits a lot more other particles than IMF, from the central-limit theorem, V th 
is expected to be close to a Gaussian function, 

r> try \ ( (Z oth — (Z oth )) 2 \ 

^(^)ocexp(^ 2a2{Zoth) y (13) 

The dispersion in Eq. (13) should be primarily associated with light (Z=l and Z=2) particles 
in the participant region, and, possibly, to some extent with the amount of charge that the 
spectator matter carries off. From the central-limit theorem, given Poisson distributions, 
the dispersion is then a 2 (Z oth )> (nz=i) + 4(nz=2)- Substituting Eq. (13) into Eq. (12) we 
get 

Pm(n)« E <E^U{P;M) exp ( JEA Zvti K)2 ) ■ ^ 

{n^-M&MF} T V K^oth) / 

Finally, to assess the effect of the constraining factor in Eq. (14) we approximate the charge 
in the exponential by its average value, z v ~ (z) = ^^(n^z^/ (n), and obtain 

P m (n) oc P p (n) exp (-a(n - (n)) 2 ) , (15) 

where a = (z) 2 /2 a 2 (Z oth ) < (z) 2 / '2((n z =i) +4(n z=2 )). Charged particle multiplicities mea- 
sured in Ref. [24] suggest a< 0.2. 

The main effect of the constraint in Eq. (15) is to narrow the multiplicity distribution, 
compared to the Poisson distribution. Depending on correlations between charge, mass, 
and energy within the rest of the system, the other constraints may affect the multiplicity 
distributions further. 

This work is supported by the National Science Foundation under Grant numbers PHY- 
95-28844 and PHY-9403666. 
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FIGURES 

FIG. 1. Three probability distributions for the constrained Poisson distributions of Eq. (8), for 
a = 0.1, A = 3 (solid points), A = 6 (open points) and A = 10 (solid points). The solid and dashed 
lines are fits with the binomial distributions, Eqs. (3-5). The corresponding fitting parameters, m 
and pb are listed in the figure. 

FIG. 2. Dependence of 1/n on l/y/a/b+(n) for a/b =0, 0.5 and 1.0. 

FIG. 3. Left panel: Arrhenius-type plot for the Ar + Au reaction at E/A = 110 MeV [2, 
4]. The solid line is the self-correlation of l/(n) as a function of 1/^/0.5 + (n) scaled according 
to Eqs. (4) and (11) with the experimental values of m=12, 6=213 MeV [2,4,24]. Right panel: 
Dependence of l/(n) on l/\fW c . 
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